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METHOD OF CALIBRATING BIAS DRIFT WITH TEMPERATURE FOR 
A VIBRATING STRUCTURE GYROSCOPE. 

This invention relates to a metliod of calibrating bias drift with operating 
temperature over an operating temperature range for a vibrating structure 
5 gyroscope. 

Vibrating structure gyroscopes are known in many forms such as tuning 
fork, cylinder or planar ring structures. Planar ring structure gyroscopes can be 
made of silicon with a silicon vibrating ring structure movably mounted by a 
number of radially compliant legs on a fixed support. Drive means are provided 
10 for putting and maintaining the ring structure in vibratory motion at the natural 
resonant mode of the stmcture. Pick-off means are included for detecting 
vibration of the vibrating structure, giving a signal proportional to velocity. 

Such a silicon planar ring structure gyroscope is relatively simple to 
manufacture in bulk with low production costs but only has a limited bias drift 

15 performance of about 1 degree per second over its typical operating 
temperature range of between -40°C and +80°C. Whilst this performance is 
satisfaQtory for civil automotive market needs it would be advantageous if the 
bias drift perfomnance could be improved to the order of between 1 to 10 
degrees per hour which would enable such an improved performance 

20 gyroscope to be utilised for other applications such as navigation. Bias drift is 
defined as the variation of the gyroscope output signal as a function of time and 
temperature in the absence of an applied rotation rate. 

A simple way of improving the performance of the bias drift over 
temperature would be to model the bias with an external temperature sensor 
25 with the use of a polynomial representation as below: 

B = a,+ a,S{T) ^a^SiTf + a^Sijf + ... - - - (1 ) 

where S(T) is the signal representing temperature from the temperature sensor 
and B the bias. In order to can^ out this calibration process the vibrating 

of time cycled over the operational temperature envelope, data-logging both the 
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bias (while the gyro is stationary) and the output of the temperature sensor. The 
rate sensor is held fixed during the period of testing. 

A multiple linear regression (such as available from typical spreadsheet 
such as Excel) can be used to determine the coefficients of the fit in a least 
5 squares fit manner from data over the operational temperature range. It is then 
assumed that the fit parameters ai, 32 etc. are constant with time, so that this 
one off calibration can then be used for the tested calibration gyroscope in 
operation. 

The problem with this temperature sensor approach is that it is not 
10 closely related to the gyroscope characteristics which accounts for the limited 
scope of improvement that can be achieved. A temperature sensor will not be 
in exactly the same physical position as the gyroscope so there will inevitably 
be temperature differences especially under dynamic temperature 
environments. This will result in either a lead or lag of the temperature sensor 
15 with respect to the gyro results in apparent hysteresis in the modelling with 
increased errors as shown in Figure 7 of the accompanying drawings, as will be 
later described 

There is thus a need for a generally improved method of calibrating bias 
drift with temperature for a vibrating structure gyroscope which at least 
20 minimises these problems. 

According to one aspect of the present invention there is provided a 
method of calibrating bias drift with operating temperature over an operating 
temperature range for a vibrating structure gyroscope having a substantially 
planar, substantially ring shaped silicon vibrating structure, primary drive means 

25 for putting and maintaining the vibrating structure in carrier mode resonance, 
and secondary drive means for nulling response mode motion of the vibrating 
structure, which secondary drive means includes means to separate a detected 
response mode motion signal into a real component induced by applied rotation 
of the vibrating structure gyroscope and a quadrature component which is an 

30 error term indicative of error mismatch between carrier mode resonance 
frequency and response mode resonance frequency, including the steps of 
measuring, over an operating temperature range of the vibrating structure 
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gyroscope, primary drive means voltage P which is a measure of change in 
quality factor Q of the vibrating structure with temperature, vibrating structure 
frequency f which is a measure of change of temperature of the vibrating 
structure, secondary drive quadrature component values Sq which is a measure 
5 of real component bias errors with temperature, and secondary drive real 
component values Sr which is a measure of change in bias, that is the zero 
inertial rate offset, of the vibrating structure gyroscope with temperature 
substituting the values obtained in the relationship 

k I m 

10 where aki„j are bias calibration coefficients for the vibrating structure gyroscope 
over the operating temperature range, and calculating from the relationship the 
coefficients auni to provide a set of bias calibration coefficients for the vibrating 
structure gyroscope over the tested operating temperature range. 

Preferably the coefficients akim are calculated from the relationship by 
15 carrying out a multiple linear regression on the relationship. 

Conveniently the coefficients a/dm are calculated by Kalman filtering. 

For a better understanding of the present invention, and to show how the 
same may be carried into effect, reference will now be made, by way of 
example, to the accompanying drawings, in which: 

20 Figure 1 is a graphic illustration of the normal mode vibration pattern for 

a silicon vibrating structure gyroscope, 

Figure 2 is a graphic illustration of the vibration pattern for a silicon 
vibrating structure gyroscope when there is misalignment between the normal 
modes and the primary drive and pick-off means by an angle oc, 

25 Figure 3 is a graphic illustration of the Corriolis Effect on the primary 

mode vibration pattern for a silicon vibrating structure gyroscope, 

Figure 4 is a diagrammatic representation of a conventional electronic 
control circuit for a conventional vibrating structure gyroscope having a silicon 
ring shaped vibrating structure. 
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Figure 5 is a block function diagram of the control circuit of Figure 4, 

Figure 6 is a block function diagram of a modified version of the control 
circuit of Figure 4, 

Figure 7 is a graphical representation of bias versus frequency showing 
5 the hysteresis effect on bias over temperature. 

Referring now to Figures 1 to 4 of the accompanying drawings a silicon 
ring shaped vibrating staicture 1 has two fundamental modes, differing by, 
typically, less than 1 Hz in frequency. Ideally the primary mode is driven at its 
natural frequency, and the secondary mode Is suppressed. In practice the ring 
10 shaped structure 1 oscillates at a frequency between the primary and 
secondary natural frequencies, and both modes are excited. The magnitude of 
oscillation is held approximately constant across the full operating temperature 
range. An AGC(automatic gain control) loop 2 operates to control the amplitude 
of the primary motion. 

15 The gyro operating principle involves the use of Coriolis forces coupling 

the primary and secondary modes when inertial rotation occurs about an axis 
nomnal to the plane of the ring shaped structure 1. The force needed to 
suppress the Coriolis force Is proportional to the inertial angular rate. 

The main error driver is the angular misalignment of the primary drive 
20 and the primary mode as shown in Figure 1. This misalignment angle varies 
across the temperature range, and is subject to hysteresis. It is not simply a 
function of present temperature, but Is a hysteretic, stochastic function of 
temperature history, which is unknown. It is this hysteresis that gives rise to 
errors that cannot be calibrated out off-line by conventional means. Apart from-.- 
25 the output of measured gyro rate 3 (corrupted by unknown errors), additional 
measurements are available from the fxilmary drive vottege 4, frequency 5, and 
secondary quadrature drive or bias 6. Thase measurements are known to be a 
function of the various error drivers. 

Figure 1 of the accompanying drawings shows the primary normal mode, 
30 with primary drive means 7 and primary pick-off means 8. Also shown, at 45 
degrees to the primary mode, are the axes of the secondary mode with 
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secondary drive means 9 and secondary pick-off means 10. The mode position 
and driver and pick-off positions are all shown in their ideal positions. The 
secondary mode is rotated 45 degrees relative to the primary. 

Misalignment of the normal modes with the drivers and pick-offs by an 
5 anglea, as shown in Figure 2, results in a combination of the primary and 
secondary modes being used to null motion at 45 degrees to the primary drive 
means 7. The primary pick-off means 8 will detect mainly the primary mode, 
but also the secondary mode, to a lesser degree. The primary drive is the 
reference direction from the gyro centre, from which all other angles are 
10 measured. 

Due to nominal symmetry of the ring shaped staicture 1 , it appears that 
the primary mode could take up any position on the structure 1, however minor 
asymmetries may determine its actual position. Thus the mode position tends to 
be locked to one angle a due to anisotropics around the ring shaped structure 
15 1. U, in Figure 2 represents the primary mode amplitude and 11 the primary 
mode node. 

Figure 3 shows the velocity vectors 2 at four points for the primary mode. 
Due to Coriolis effect in the presence of rotation, there are D'Alembert's forces 
13 applied at each point. Cleariy, the primary mode motion in the presence of 
20 rotation gives rise to forces, which drive the secondary mode, and by symmetry, 
the secondary mode motion in the presence of rotation gives rise to forces, 
which drive the primary mode. 

The pattern of Coriolis acceleration suggests that the secondary mode 
will be driven by D'Alembert's forces. This leads to a cross-coupling between 
25 the primary and secondary normal modes, in both directions. There is a gain 
associated with this effect, the Bryan gain. Thus, for a rotating ring structure 
there is a special effect, referred to here as the Bryan Effect, which forms the 
physical basis of the gyro's operation. 

Typically a silicon ring shaped vibrating structure gyroscope is controlled 
30 by an electronic control system as shown in Figure 4 of the accompanying 
drawings. This has a primary loop 14 and a secondary loop 15. Making up the 
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. primary loop 14 is a frequency control loop 16 and an amplitude control loop 17. 
The loop 16 contains a phase locked loop 18 and a voltage controlled oscillator 
19. The loop 17 contains an automatic gain control 20 to stabilise the amplitude 
of motion, having an input V(agc) voltage to set the amplitude of motion. Loop 14 
5 is connected between the primary picl<-off means 8 and the primary drive 
means 7. The secondary loop 15 is connected between the secondary pick-off 
means 10 and the secondary drive means 9 and acts as a 
demodulating/remodulating null control loop. To this end it contains a 
quadrature component loop 20 and a real component loop 21. 

10 Figure 5 shows a block function diagram of the control system of Figure 

4, in which like parts have been given like reference numbers. The output from 
the automatic gain control 20 is indicated at 23. 

Figure 6 shows a block function diagram of the control system of Figure 
4, in which like parts have been given like reference numbers, but also in which 
15 there is an extra servo loop 24 connecting the phase locked loop 18 to the 
gyroscope vibrating structure 1 to ensure that there is a 90'' phase shift between 
the primary drive voltage 4 and the primary pick-off signal. 

Over temperature, due to the change of modulus of the silicon and the 
size of the ring structure 1, there is a change of resonance frequency of 

20 -0.4HZ/C. With a typical phase error of the loop of 0.5 degrees, the temperature 
of the ring structure can be determined to an accuracy of 0.01Hz corresponding 
to 0.025 C. The frequency of vibration of the structure 1 can be determined 
very accurately during calibration and this implies that the temperature 
environment of the ring structure can be very precisely determined with no time 

25 lag between the "temperature sensor" (the frequency of the ring) and the ring 
structure itself. 

The frequency of the ring structure 1 when It is vibrating at resonance is 
a good measure of the actual temperature of the ring structure with a typical 
coefficient of -0.4Hz/C. The frequency and the gyro bias with the gyroscope 
30 stationary in a temperature chamber can be recorded over the operational 
range, which is typically — 40C and an improved model for the gyro bias can be 
written as: 
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where f is the frequency of the ring. This again is a calibration which is carried 
out over the operational temperature range of the gyro with data being collected 
for the gyro bias b, scale factor s and frequency f of the ring structure as a 
function of time when a typical temperature profile as in Figure 7 is found. 

Whilst this model results in an improvement in the gyro bias performance 
over the external temperature recording process, it cannot eliminate the effects 
of hysteresis of bias over temperature as there is no time history in the model. 

Figure 7 shows that for variation in temperature there is a hysteresis 
effect in that the solid line in the graph shows a deviation depending upon 
whether the temperature is increasing or decreasing, as shown by the arrows, 
from a mean track 25. Hence this process for calibrating bias drift is inherently 
inaccurate. 

With reference to the block diagram in Figure 6, there are a number of 
measurands in the gyro as follows: 

a) primary drive voltage 4, given the symbol P 

b) Frequency 5, given the symbol f 

c) Secondary quadrature drive levels 6, given the symbol Sq, and 

d) Secondary real drive 3, which is the normal rate output Sr 

The primary drive voltage P is a measure of the drive required to maintain the 
primary oscillation of vibrating structure 1 at a fixed amplitude. As the 
temperature changes the Q of the gyroscope decreases with increasing 
temperature requiring a higher drive level. Thus the primary drive voltage P is a 
measure of 1/Q as a function of temperature, and this knowledge can be used 
to take out the effects of temperature variation of Q. The frequency, as 
described before, is a measure of ring structure temperature and represents a 
very accurate indicator. The secondary quadrature signal Sq is related to error 
effects of the real bias. In particular the ring shaped vibrating structure 1 is 
known to have two eigen modes . 
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These modes are not naturally aligned to the primary and secondary 
drive axes, but have a misalignment angle a. When the primary drive is 
energised, if a is not equal to zero then a linear combination of the two eigen 
modes is excited. These eigen modes are characterised by resonant 
frequencies fi and f2 and quality factors Qi and Q2. If the two modes do not 
have the same resonant frequency, and the angle a is not equal to zero then a 
quadrature drive signal is required to null out the quadrature motion. It can be 
shown that the quadrature level is related to the mode angle a. 

There are also errors in the bias , the zero inertial rate offset of the gyro, 
which is measured by Sr. it is known that when the Q factors of the two modes 
are different, then there is a bias error proportional to sin 4a and the values of 
AQ. It is seen that there is thus a correlation between the value of Sq and the 
proportion of the error in Sr due to AQ. Thus the variation of Sr will track the 
variation of Sq, so that the error due to delta Q in gyro bias can be 
compensated. 

Over the temperature range of ±80° this may be written as: 



where an is a fit coefficient from a least squares fit. 

In a similar way the gyro bias will have error terms which vary as the Q of the 
mode changes, which is measured by the primary drive level, and this can be 



... (4) 



written: 



---(5) 
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Putting all these terms together we can write: 

'5.=Z/"S^;S^'X. "-(6) 

Note that in this expression the values of k, I and m can be both positive 
and negative. 

5 In order to cover the major errors that are known on the gyro typically a 

total of 7 terms are used . 

These cover errors such as: 

a) phase error of the primary drive with respect to the secondary pick off. 
These are electronic phase errors with a typical value of 0.5 degrees due 

10 to limited pre arnplifier bandwidth and time delays through integrated 

circuits. 

b) Effects of delta Q as discussed above 

c) Effects of the angle offset between the secondary pick off and primary 
drive 

15 d) Effects of bias voltage to null the original bias offset 

e) Effects of inductive pick up between primary drive and secondary pick off 

f) Effects of capacitive pick up between primary drive and secondary pick 
off 

g) Effects of magnetic field anomalies. 

20 By a suitable choice of the powers of k, I and m it is possible to model 

these 7 error terms using the three measurands. During the calibration process 
the values of the measurands - Sr Sq P and f are data-logged as the 
temperature varies over a range of ± 80''C. A multiple linear regression is 
carried out to determine the coefficients akim- 

25 Other analysis techniques can be used to derive the coefficients from the 

data such as Kalman filtering techniques. 
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As the measurands are correlated to the gyro bias errors drivers, 
measurements of the these terms will improve the bias compensation. The 
polynomial form given above is then linked back to physical error sources (as 
discussed), and the closer the calibration is to the physical reality then the 
5 better the modelling that can be achieved. 

As a result It is now is possible to model gyro bias hysteresis, which a 
simple polynomial model cannot achieve. This is due to the correlation of real 
and quadrature bias in particular. Therefore if there is hysteresis in real bias 
there will generally be a corresponding hysteresis in quad bias as the error 
10 mechanisms (such as variation of alpha and hysteresis of alpha with 
temperature) are common to these two. 

In practice the calibration coefficients awm are established for a specific 
gyroscope by applying in an oven, a desired temperature sequence to the 
gyroscope. Test data, including the previously specified measurands is 
15 recorded and a set of bias calibration coefficients established for the gyroscope. 



